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Using first-principles density functional theory calculations, we discover an anomalously large 
bi-axial strain-induced octahedral rotation axis reorientation in orthorhombic perovskites with 
tendency towards rhombohedral symmetry. The transition between crystallographically equivalent 
(isosymmetric) structures with different octahedral rotation magnitudes originates from strong strain- 
octahedral rotation coupling available to perovskites and the energetic hierarchy among competing 
octahedral tilt patterns. By elucidating these criteria, we suggest many functional perovskites would 
exhibit the transition in thin film form, thus offering a new landscape in which to tailor highly 
anisotropic electronic responses. 



PACS numbers: 61.50.Ks, 68.55.-a, 77.55.Px 

Introduction. Phase transitions are ubiquitous in na- 
ture; they describe diverse topics ranging from crystal- 
lization and growth to superconducting Cooper pair con- 
densation. Isosymmetric phase transitions (IPT) — those 
which show no change in occupied Wyckoff positions or 
crystallographic space group — are an intriguing class since 
there are relatively few examples in crystalline matter 
PQ: Most condensed matter systems respond to external 
pressures and temperatures by undergoing "conventional" 
symmetry- lowering displacive [2j , martensitic [3] or recon- 
structive [1] transitions. Furthermore, the experimental 
characterization and identification of a suitable symmetry- 
preserving order parameter through such transitions is 
often challenging [5]. Although some electronic order pa- 
rameters [5] [7] that include ferroelectric [BHTT] or orbital 
polarizations [H] have been proposed for IPT, which lead 
to subsequent changes in local cation coordinations [TBr - 
IT5] , to the best of our knowledge, there is no case where 
the IPT connects two structures with essentially the same 
local bonding environment. 

Using first-principles density functional calculations, 
we find an isosymmetric transition in the low energy 
rhombohedral phases of epitaxially strained orthorhom- 
bic perovskites and describe how to experimentally ac- 
cess it. We show that the transition originates from 
non-polar distortions that describe the geometric connec- 
tivity and relative phase of the BOq octahedra found 
in perovskites. Although a previous IPT in a thin film 
perovskite that relies on strong strain-polar phonon cou- 
pling has been reported [16], we describe here a universal 
symmetry preserving transition that originates from the 
strong lattice-octahedral rotation coupling ubiquitous in 
nearly all perovskites. For this reason, the large isosym- 
metric reorientation of the oxygen rotation axes should be 
readily observable in many rhombohedral perovskites with 
diverse chemistries. Since the dielectric anisotropy in per- 



ovskites is strongly linked to deviations in the octahedral 
rotation axis direction [17j . we suggest control over this 
transition could provide for highly tunable high-K dielec- 
tric actuators and temperature-free relative permittivity 
resonance frequencies |18j . 

We choose LaGa03 as our model system since it has 
a tolerance factor of r = 0.966 indicating the perovskite 
structure is highly susceptible to GaC>6 octahedral rota- 
tions about the principle symmetry axes 19J: At room 
temperature it is orthorhombic Pbnm and undergoes a 
first-order phase transition to rhombohedral i?3c around 
418 K [20], with a subsequent change in the GaOe octahe- 
dral rotation patterns from a~a~c + to a~a~a~ , respec- 
tively, in Glazer notation [21]. The + (— ) superscripts 
indicate in- (out-of)-phase rotations of adjacent octahe- 
dra along a given Cartesian direction. The non-magnetic 
Ga 3+ cations additionally allow us to eliminate possible 
contributions of spin and orbital degrees of freedom for 
driving the IPT through electronic mechanisms. 

Calculation details & notation. Our density functional 
calculations are performed within the local density ap- 
proximation (LDA) as implemented in the Vienna Ab 
initio Simulation Package (vasp) [22] 123 with the pro- 
jector augmented wave (PAW) method [24], a 5 x 5 x 5 
Monkhorst-Pack fc-point mesh [25] and a 500 eV plane 
wave cutoff. We relax the atomic positions (forces to be 
less than 0.1 meV A -1 ) and the out-of-plane c-axis lattice 
constants for the strained films [26 . 

The principle difference between the ground state or- 
thorhombic Pbnm and metastable [12 meV per formula 
unit (f.u.) higher in energy] rhombohedral R3c phases 
of LaGaOa is that the GaOe octahedra rotate in-phase 
(+) along the Cartesian z-direction of the Pbnm struc- 
ture while they rotate out-of-phase (— ) about that same 
direction in the i?3c structure. Our homoepitaxial bi- 
axial strain calculations simulate film growth on a cubic 
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FIG. 1. (Color) Evolution of the total energy (a) for the e- 
Pbnm and C'2/c phases with in- and out-of-phase octahedral 
rotations (inset) along the z-direction. (b) Magnified region 
about the IPT (shaded), (c) The change in axial ratio with 
strain shows a discontinuity in the C2/c phase that is absent 
in the e-Pbnm structure. 



(OOl)-terminated substrate. We therefore choose the c + 
rotations of the orthorhombic phase to be about the axis 
perpendicular to the epitaxial plane (Fig. [T]), to evaluate 
the bi-axial strain effect on the in- versus out-of-phase 
GaO*6 rotations present in the two phases. Note, the bi- 
axial constraint preserves the orthorhombic symmetry in 
the a~a~c + phase, however, we designate the epitaxially 
(e) strained phase as e-Pbnm to distinguish it from the 
bulk structure. In contrast, the symmetry of the bulk 
rhombohedral phase is lowered to monoclinic (space group 
C2 1 c) and we therefore refer to it as such [57] . 

Strain- stabilized structures. We first compute the evo- 
lution in the total energy with bi-axial strain for the 
e-Pbnm and C2/c structures [Fig. [it a)]. We find that be- 
tween approximately —1 to +3% strain, the orthorhombic 
phase with the a~a~c + rotation pattern is more stable 
than the monoclinic oTarc" structure. For now we fo- 
cus on the monoclinic phases [Fig. |ljb)] near 0% strain: 
We find an abrupt discontinuity in the first derivative of 
the total energy with strain for the monoclinic structure 
between two states with the same symmetry, denoted 
C2/c(l) and C2/c(2). In contrast, we find a single contin- 
uous equation of state with uniform hydrostatic pressure 
(over ±50 GPa) for the bulk structures. The evolution in 
the c/a axial ratio for these structures is also qualitatively 
different [Fig. [l|c)] . The e-Pbnm axial ratio continuously 
decreases with increasing tensile strain (consistent with 
elastic theory), whereas in the C2/c structures a sharp 
discontinuity occurs in the vicinity of c/a ~ 1. We find 
the first-order phase transition occurs at a critical strain 
of ~ 0.18% from intersection of the quadratic fits to the 
total energies. 

Microscopic structure evolution. To investigate if the 
C2/c(l) — > C2/c(2) transition is indeed isosymmetric, we 
evaluate how the internal structural parameters - octha- 



hedral tilts and bond distortions - evolve with epitaxial 
strain [Fig. |2ja)]. We find a continuous evolution in the 
GaC>6 rotation angles for the e-Pbnm structures (open 
symbols): the rotation axis changes from being along the 
[001]-direction to mainly in-plane along [110] as the strain 
state changes from compressive to tensile. In contrast, we 
find an abrupt change in the octahedral rotation angles 
with strain in the monoclinic phases (filled symbols). We 
identify that the C2/c(l) and C2/c(2) phases, despite pos- 
sessing the same symmetry are distinguishable — they have 
either mainly [110] in-plane or [001] out-of- plane GaOg 
octahedra rotations. Consistent with the orthorhombic 
case we find that increasing tensile strain drives the the 
octahedral rotation axis into the [110]-epitaxial plane. 

The bi-axial strain is not solely accommodated by rigid 
octahedral rotations. It produces additional deviations 
in the Ga-0 bond lengths and causes La cation dis- 
placements. We quantify the former effect through the 
octahedral distortion parameter A = g J2 n =i 6t(^( n ) ~~ 
(6))/ (S)} 2 , where d is a Ga-0 bond length and (5) is the 
mean bond length in the GaOg octahedra. With increas- 
ing strain, A increases, indicating that bond stretching 
(and compression) occurs simultaneously with changes 
in the magnitude of the octahedral rotation angles to 
alleviate the substrate-induced strain [Fig. |2fb)]. Ac- 
cording to our bond-valence calculations, the Ga-0 bond 
stretching alleviates the "chemical strain" imposed on 
the over-bonded Ga 3+ cations when a regular GaOe oc- 
tahedra (A — > 0) occurs. The IPT allows the monoclinic 
phase to maintain a uniform charge density distribution 
with the a~a~c~ tilt pattern; this is assisted by the anti- 
parallel La displacements [Fig. |2][c)], which change sign 
across the transition (shaded), maintaining a trigonal pla- 
nar configuration in the GslOq rotation-created cavities. 
Note, this chemically over-bonded structure is absent in 
the e-Pbnm structure because the a~a~c + tilt pattern 
(D2h symmetry) permits non-uniform Ga-0 bonds. Since 
the rotation pattern never reverses, a single La cation 
displacement direction occurs. 

Origin of the IPT. To identify the origin of the isosym- 
metric transition, we first analyze the energy of the 
monoclinic a~a~c~ structures under different bi-axial 
strain states as a function of direction and magnitude 
of the GaOg octahedron rotation axis. The direction 
of the GaOg rotation axis with the a~a~c~ pattern is 
constrained to be in the (TlO)-plane because the rota- 
tion pattern can be decomposed into a~a~c° and a°a°c~ 
rotations with with axes aligned along the [110]- and 
[001]-directions [Fig. [5Ja)]. We show in Fig. [3^a-c) our 
first-principles results of the energy dependence on the 
direction (vertical axes) and magnitude (horizontal axes) 
of the GaOg rotation axis for strain values of -1.5%, 0.0% 
and 1.5%, respectively. For all strain states, we find a 
single well-defined energy minimum for each direction of 
the GaOg rotation axis [Fig. [3ja-c)]. We therefore are 
able to remove the rotation angle magnitude as a variable 



FIG. 2. (Color) Evolution in (a) the GaC>6 rotation angles about different directions relative to the substrate, (b) the octahedral 
distortion parameter A, and (c) the La displacements about the bulk structures with epitaxial strain. 



and to analyze the energy dependence solely in terms of 
the bi-axial strain and the GaOg rotation axis direction 

We show in Fig. [3]jd) the calculated evolution of the 
extremal octahedra rotation axis directions with bi-axial 
strain: local energy minima (maxima) are indicated with a 
heavy (broken) line, and for -1.5%, 0.0% and 1.5% bi-axial 
strains, we explicitly mark the extrema using symbols. 
Consistent with our earlier structural analysis, we find 
that the rotation axis direction smoothly approaches the 
[110]- ([001]-) direction for large tensile (compressive) 
strains. For the range of strains between —0.5% and 
0.5%, we observe the co-existence of two energy minima 
separated by an energy maximum (broken line); this 
indicates an inaccessible region of rotation axis directions 
close to [111] for any value of strain and is consistent with 
a first-order transition. 

Our results suggest there are two main reasons for the 
appearance of the isosymmetric transition in epitaxially 
strained rhombohedral perovskites. The first reason is 
that the octahedral rotations are strongly coupled to the 
bi-axial strain. This coupling originates from the rigidity 
of the GaC-6 octahedra, since the rigidity causes contrac- 
tion of the crystal lattice in the direction orthogonal to 
the rotation axis [5^1 ED] . Second, the bulk rhombohedral 
a~a~a~ structure of LaGaOa is higher in energy than 
the bulk orthorhombic a~ a~ c + structure. 

We now show that the energy ordering of the bulk 
phases is responsible for the inaccessible region of rota- 
tion axis directions. The a~a~a~ and a~a~c + structures 
differ only in the phase of the GaOe octahedra rotations 
about the z-axis. Thus, each structure can be transformed 
into the other through a combination of rigid octahedral 
distortions. One distortion should deactivate the a~ ro- 
tation about the z-axis, while the other would induce 
the c + rotation about the same axis. We would expect 
these distortions to impose minor energetic penalties since 
they are nearly rigid [31) . In the present case, where the 
a~a~a~ structure is higher in energy than a~a~c + , we 
expect that introduction of either of these distortions into 
the higher energy a~a~a~ structure will lower the total 
energy [3"2] . 

Finally, smoothness of the total energy as a function 
of strain and rotation axis direction requires that the 



difference between the number of energy minima (iV m j n ) 
and maxima (N max ), for any value of bi-axial strain, re- 
mains fixed as elaborated in Morse theory j^Sj. In other 
words, any smooth deformation which produces addi- 
tional energy maximum must also produce additional 
energy minimum. Because of the first reason for the IPT 
mentioned above, we anticipate that for sufficiently large 
compressive or tensile strains, the strain-octahedral rota- 
tion direction coupling dominates to yield a single energy 
minimum: N m { n = 1, N max = and then from continuity, 
Nm\ n — N max = 1 must remain constant for all strains. 
From our energetic hierarchy of the bulk structures, we 
conclude that when strain induces structural distortions 
with magnitudes which nearly coincide with those of the 
bulk a~a~a~ phase (near 0% strain and [111] direction), 
there will exist an energy maximum (N max = 1); from con- 
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FIG. 3. (Color) Calculated energy (a-c) of the monoclinic 
LaGa03 phases as a function of the Ga06 rotation axis direc- 
tion and angle magnitude at -1.5%, 0.0% and 1.5% strains, 
respectively, with contours at 5 and 0.5 meV/f.u. (d) Position 
of energy minimum (solid line) or maximum (broken line) 
with strain; circles correspond to minima (a-c), and the cross 
indicates the saddle point in (b). 
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tinuity, this must introduce two energy minima (iV ro i n = 2) 
at the same value of strain. These reasons together pro- 
duce the energy landscape shown in Fig. |3^d) and require 
an IPT in the LaGaC>3 system. For comparison, the or- 
thorhombic phase of LaGaC>3 does not show an IPT as 
one varies bi-axial strain, since the second condition for 
the transition described above does not apply, ie. the 
orthorhombic structure is the global ground state: N m - ln 
is fixed to 1 (N max = 0) for all strain values. 

Accessing & applications of the IPT. We obtain a 
C2/c — > e-Pbmn transition near —1% compressive strain 
with respect to our hypothetical cubic LaGaC>3 phase. 
In the vicinity of the IPT at K, however, the e-Pbmn 
phase is the global ground state. Although our minimal 
model for the IPT relies on this energetic ordering of the 
competing rotational phases (a~a~a~ versus a~a~c + ), 
we anticipate three experimental routes by which to access 
the essential signature of the isosymmetric transition — 
large strain-induced reorientation of the octahedral rota- 
tion axis direction. First, the monoclinic phases could be 
stabilized in thin films through the substrate coherency 
effect [33] , where the film's tilt pattern adopts that of 
the substrate: Perovskite substrates with the a~a~a~ 
(LaAlOa) or the a°a°c~ (tetragonal-SrTiOa) tilt pattern 
are promising candidates. Second, additional electronic 
degrees of freedom (first- and second-order Jahn- Teller 
effects), introduced through cation substitution, could be 
exploited to stabilize the IPT because they often energet- 
ically compete with the octahedra rotations [34]. Lastly, 
experiments performed above the bulk LaGaC>3 structural 
transition temperature (~ 100°C) would make the mon- 
oclinic phases accessible at all strains. The IPT would 
exhibit a weak-first order transition while still providing 
strong strain-octahedral rotation axis direction coupling. 
At sufficiently high temperatures, the IPT could be sup- 
pressed and its boundary terminated by a critical point 
[35]. 

We have shown that strain-octahedral rotation axis 
directions are strongly coupled in epitaxial perovskite thin 
films. We suggest similar large reorientations of coordinat- 
ing polyhedra frameworks could be achieved in alternative 
structural families: thin films with the garnet, apatite 
or spinel structures are particularly promising. However, 
the functional materials design challenge remains: how 
does one couple the rotation axis direction to additional 
electronic degrees of freedom? For this reason, we advo- 
cate for detailed epitaxial film studies on perovskites close 
to the i?3c -f-> Pnma phase transition (0.96 < r < 1.01). 
Controlling the IPT in LaCr03, LaNi03 and LaCu03 per- 
ovskites could yield unknown, and potentially functional, 
orbitally-, spin- and charged-ordered phases. 
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